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Abstract 

When rotation is not taken into account, the measurement of the Gravitational 
Redshift can provide unique information about the compactness {M/R) of the star. 
Rotation alters the gravitational redshift rendering thereby a unique determination of 
the compactness parameter impossible. Nevertheless, it can be shown that by using 
some theoretical input, useful information on, say, the radii of compact rotating 
objects can still be extracted. Moreover, by measuring the gravitational redshift 
one can infer the maximum angular velocity of the object. As it is well known, the 
minimum observed periods of rotation are found in millisecond pulsars. Here we 
show that millisecond periods are actually a semi-theoretical limit that can be found 
by General Relativistic arguments corresponding to the maximum angular velocity. 
We apply our method to compact objects such as pulsars, white dwarfs and neutron 
stars. 

Keywords: maximal angular velocity; stars: neutron; stars: rotation 
PACS: 04.40.Dg, 97.60.Jd 

1 Introduction 

Compact Objects such as neutron stars [1] [5] have radii very close to their Schwarzschild 
radii and hence General Relativity should be used to describe the gravity close to their 
surface. Neglecting rotation, the geometry of space-time can be described using the well 
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known spherically symmetric Schwarzschild geometry and information on the ratio M/ R 
of a compact object can be obtained from the (observed) gravitational redshift ([3], [4]). 
This is a useful, model-independent way, to gain insight into properties of neutron stars 
and white dwarfs. Although it is a model-independent procedure, it relies on certain 
basic assumptions, like the absence of rotation and a perfect spherical symmetry (i.e. 
zero quadrupole moment). To gain a deeper understanding it makes then sense to relax 
some of these assumptions and to examine what information can be extracted from the 
gravitational redshift in the more general case. In the present paper, we shall study 
the situation with a non-zero angular velocity SI, but assuming that the object is rigid 
enough to allow the approximation of a negligible quadrupole moment. In particular, 
when rotation is taken into account, spherical symmetry is lost and off diagonal terms 
appear in the metric which has the following general form: 

ds^ = gttdf + grrdr^ + geedO"^ + g^.^.dcj)'^ + 2g^td(j)dt. (1) 

The exact form of this metric including magnetic field, quadrupole moment and even 
radiation ([5], [5], [7]) is still a subject of research. We will use a first and second order 
approximation to this metric neglecting the effects of the quadrupole moment and the 
magnetic field. Furthermore, it is expected that only relatively young Neutron Stars can 
have a differential rotation [T] (encountered also in boson stars [S]) which we therefore 
neglect here. 

Deriving the gravitational redshift from ^ we will show that, for a fixed mass, a 
unique solution for the radius of the compact object does not exist. Instead, we obtain 
two different solutions which can differ by orders of magnitude as long as the angular 
velocity is not too large. In this case, all theoretical models still favor the solution of 
the smaller radius and we can select this solution without the necessity to refer to some 
details of a specific model. With increasing angular momentum, the two radii approach 
each other and the above selection rule is not effective anymore. However, in such a case 
we can define a narrow range of allowed radii which is still a valuable model independent 
information for fast rotating objects. 

As a bonus of our examination of the properties of the gravitational redshift Z, we can 
show that there exists an upper bound on the angular velocity depending on the redshift. 
One can even put an absolute upper bound which is independent of Z. Interestingly the 
bounds come out to be in the range of the observed millisecond pulsars. This in turn 
allows us to assume that the angular velocity of millisecond pulsars is indeed approaching 
its maximal value. The consequences of this assumption will be discussed in the text 
below. One of them is the confirmation of the result that the light emitted in very fast 
millisecond pulsars stems mostly from the equatorial region ([9], [10], [11]). We compare 
our results to the so called 'Mass Shedding Limit' [T2] and to other work related to the 
maximum angular velocity. 

One of the main results of the paper, discussed in section four, is that our analytical 
formula on maximal angular velocity is comparable to results obtained by using extensive 
numerical calculations. In particular, we agree with the numerical findings of reference 
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The essential parameters which enter our equations are the radius R, the angular 
velocity and the mass M . The minimum central density at which a neutron star is stable 
is simply the density at which neutrons become unstable to beta decay (po ~ 8 x 10^ g/cm^ 
[13]). Using the well known Oppenheimer-Volkov [T^ equations and a plausible equation 
of state, one can construct a stellar model which provides the minimum mass of the 
neutron star to be around O.O8M0 which is a bit unrealistic taking into account that 
neutron stars are remnants of supernova explosions [16j . A more realistic value for the 
minimum mass is of the order of « IM© [T7], which is actually closer the maximum 

mass of a white dwarf (1.44M0). The maximum mass of a neutron star can be found 
from causality arguments [18], by recalling that the speed of sound in dense matter has 
to be less than the speed of light < c^). This condition gives a maximum mass 
of -M™s'^ ~ 3M0. The radii corresponding to the maximum and minimum masses can 
be found using the Oppenheimer-Volkov equation. Assuming an equation of state for 
a degenerate neutron fluid, the corresponding radii for the two extreme masses lie in 
the range ^ 10 km to ~ lOOfcm. This radius range can of course change if one uses a 
more sophisticated equation of state, subject of current controversy [17], but the orders of 
magnitude remain the same [19j . The canonical neutron star mass and radius are thought 
to be ~ IM© and ^ 10 km. A useful observational quantity, which agrees in order of 
magnitude with theoretical predictions, is the mean "measured" mass of the neutron stars 
in a Gaussian ensemble [20j . namely 



This value will be used if no other information on the masses is available. 

The paper is organized as follows. In the second section we will briefly discuss the 
gravitational redshift as emerging from a perturbative axial symmetric metric. In section 
three we will use these results to determine the radius or the range of the radii of the 
compact object. Section four is devoted to the maximal angular velocity derived within 
our approach. In section five we discuss some improvements by taking into account more 
terms in the expansion of the metric. In section six we apply our results to some chosen 
compact objects like a white dwarf and neutron stars. Finally, in section seven we present 
our conclusions. 

2 The Gravitational Redshift in a Perturbation Ap- 
proach 

A far away observer can measure a pulsar's angular velocity Q, given by 

dt dr dt u*' ^ ' 

Using ([3]), the four velocity of a stationary point on the surface can be written as 



(A/)ns = (1.35±0.05)M( 



(2) 




(4) 
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Through the normahzation condition of the four velocity {u^u^^ = —1), we obtain the 
time-hke component of the four velocity in terms of the metric ([ij and angular velocity 

u* = (-.9tt - 25t0f^ - ^'^g4>4>) ■ (5) 

The redshift factor can be calculated simply by recalling that the energy of radial a photon 
(?i = c = 1) is simply [3T] 

,,dxu f [ dt d(l)\ 1 ( dt d6\ , , 

with A an affine parameter and E and L conserved quantities due to the existence of two 
killing vectors. As a consequence we can write for lu 

E + QL 

^ = 172 ' (^) 

{-gtt - '2gt<p^ - ^'^9(i><p) 

and the energy perceived by a distant observer uj can be now expressed through 

ll> = Zll>o, (8) 

were uiq is the energy at the surface and Z is the redshift factor. Explicitly, the latter is 
given by [22] 

Z = {~gu - 2gt4,n - fl^g^^) . (9) 

This redshift can actually be measured for many objects [24]. Indeed, at the end of the 
paper we will employ the results of such observations. The behavior of Z as a function 
of the radius R for a fixed mass is displayed in figure la. 

The calculation of Z can be made more concrete when we consider the metric ([T]). 
To a first approximation, Zeldovich and Novikov [27] apply small perturbations to the 
Schwarzschild geometry. The elements of the metric ([T]) can be calculated to be (27j 

/ 2GM\ 2GJsiT?e R^sin^e 

5« = - (^1 - , ff*^ = — R^' = 

Here, the condition for slow rotation is given by J ^ MRgC {Rg being the Schwarzschild 
radius) [27]. Equation ([TO)l allows us to calculate Z explicitly. One obtains 

2(M,i?,17,^).(l--^ -R ^) ' (11) 

Throughout the paper we will be using the Newtonian approximation for the angular 
momentum, (J = |Mi?^ri) [28]. In view of the results obtained in [29], this is a well-based 
assumption violated only for extremely high angular velocities (which are not exceeded 
here). Taking this into account, equation (jlip simplifies to 
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Z (M, R,Vl!)=\\- a-^ - [3 MR n'^ - 7 R'^fl''^ ] , (12) 
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Figure 1: The top figure 'a', shows the behavior of the redshift equations as a function 
of the radius for different angular velocities. The bottom figure 'b' shows the behavior of 
the corresponding polynomial pS]) in units of meters. In both figures M = 1.44M0. Note 
how the gap between the two solutions becomes narrow with increasing angular velocity. 



were a, (3 and 7 are constants given by 

2G ^ 8G 1 
"=^'/^=^'7 = ^ (13) 

and we have absorbed sin 9 into the angular momentum by defining 

n' = smen (14) 

The equations ^ [T31 [H]) can be now used to either solve them for the radius by assuming 
a mean mass or a mass range and a measured angular velocity or, alternatively to predict 
the redshift. Both ways will be used below. 

3 Determination of the Radius 

It is possible to take two different approaches when using equation (fT2|) . The first one is 
to demand that the term inside the parenthesis of (fT2|) should be greater than zero, such 
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that after factoring out 1/\/R one arrives at 

/(O; R, M, Vt) = R~aM - (i MR^ 17^-7 R^^f > 0. (15) 

The hmiting values of R correspond to the equal sign in the above equation. Since this 
does not depend on Z, these values have an absolute character in the sense that they 
give the maximal and minimal radius for any compact object with mass M and angular 
velocity f2 regardless of the value of Z. Similar reasoning applies to any other quantity 
derived from (jlSp (e.g. ri[-,jax ™ ^^'^ next section). 

The behavior of the function /(O; R, M, fl) versus R is shown in figure lb. The figure 
displays the global properties of this function (which can be also inferred easily analyti- 
cally), like the local maximum and the two zeros, one of them close to the Schwarzschild 
radius. 



On the other hand, we can solve the following cubic equation for the radii 

/(Z; R, M, n) = {1- Z^)R -aM - P MR^ 17^ _ ^ ^3^2 ^ q (-j^g-j 

Obviously, this is is the same equation as (fT5|) if we put Z to zero in (fT6|) . Hence, we 
can continue examining equation ()16p and discuss the absolute limits by putting Z = Q 
at the end. The function / with non-zero gravitational redshift has the same global 
properties as The solutions of (fTB)) can be obtained analytically by parameterizing 

the Cardano formulae [30] . By a simple transformation one can get rid of the quadratic 
term in the cubic equation arriving at + py + q = 0. Depending on the sign of the 
discriminant D = (p/3)'^ -I- {q/2)'^, one can parametrize the solution using the auxiliary 
variable T — sgn (q) a/|p|/3. In our case D > and we parametrize the solutions through 
the angle a given by cos a = qjeT^ The analytical solutions are then 



were 



and 



n = -A-^cos<{ y + icos ixj' (17) 
ra = -A-^cos<{ y + icos-ix^ (18) 



A-^, ^^^\/T + r^. (19) 
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27 I 1024 (GM^X^ j_ 9GMO ^ 8(1-Z^) GMil \ 



64 ( GMQ. ^ 
25 V I 



3(1-Z2)} 



3/2 



(20) 



In figure 2b we have plotted r\ and r2 versus the angular velocity f2. The upper branch 
of each curve corresponds to the bigger radius meeting the lower value at some £7 (this 
will be discussed in the next section in more detail). With growing the difference Ar 
between the two solutions becomes smaller, however, at relative small angular velocity 
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Figure 2: In these figures we show the behavior of the two solutions ri and r2 as a function 
of the angular velocity for different masses. The top figure 'a' corresponds to the numerical 
solutions using the extended metric ([55]) . while the bottom figure 'b' corresponds to the 
analytical solutions P?! [T5|) . Note how the solutions meet at a particular angular velocity 
for each mass. 
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we can opt safely for the lower value of the obtained radius (close to the Schwarzschild 
radius) as this is favored by models. Such a point of view is not possible anymore with 
increasing fi. Note also that the solutions ri_2 do not determine a "range" for the radius 
in the strict sense. However, since both solutions depend on the mass, a range in the mass 
will determine a range for each of the two solutions. With this in mind, we can define a 
narrow range for the radius which will be explained below. 

A plausible range for the mass can be given by the close extremes in the Gaussian 
distribution ^ as discussed in the introduction. For every limiting mass we will get a 
curve like in figure 2, say Ci for the lower mass and C2 for the upper value. The curves 
with an in-between mass will fill the space between the curves Ci and C2. A vertical 
tangent to the cusp [31] of Ci will intersect the curve C2 in two points, roi and ro2, which 
we can take as a definition of range of allowed radii. The result is a single narrow range 
which becomes smaller with increasing angular velocity and is zero at a maximum fl. For 
instance, in the case of the neutron star PSR B1937-I-21, roi — 17600 m and ro2 = 25800 
such that Ar = 8200 m. The existence of the maximum angular velocity corresponding 
to the cups of every curve allows even to sharpen this concept to be discussed in the next 
section. 

It is of some interest to expand these solutions (ri and 7^2 ) neglecting small terms [52] 
in A and x- The relevant quantities can now be approximated as 



, 8GM ^ 2c /I „ 



and 



^ C3(l- ^2)3/2 5^3^13^ 

When X ^ 0, one can express ri and r2 as 



(22) 



n » ~A^A\^~^~ (23) 
\2 6 12V3 / 

r2 « -A + a(^ + ^], (24) 
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The final results of our approximation reads 
ri ; 



1 / 1 y2 _ GM _ 12G'M _ G^M^njAZ^ - 19)^ 

~ ~ C2(l-Z2) " ~15? 150c5(l-Z2)5/2 ^ > 



2GM 8G3m3172(i9_ 4^2)3 

c2(l-Z2)+ 3375c8(l-Z2)4 

Note that the first term in r2 is the Schwarzschild radius modified by a factor of 
l/(c2(l — Z2)). Actually, this term is the same result one would obtain by using a Spher- 
ically symmetric metric with no rotation. 

It is important to remember that this approximation starts to fail when % — )■ 1, which 
occurs at Z 0. 



8 



4 The Limiting Angular Velocity 



To understand the origin of a maximal angular velocity it is instructive to look at the 
generic behavior of the function /. As already briefly mentioned, the latter will have two 
zeros on the positive axis and a local maximum between them. Obviously, the case where 
the local maximum falls below zero is a limiting case corresponding mathematically to 
D = or alternatively to ri — r2 (or in a yet different method setting x = 1 (eq. 1^0)) and 
physically corresponding to a maximally allowed angular velocity). The cusps |31| of the 
curves in figures (2a and 2b) display this behavior. After some algebra one obtains 

^'maAZ) = - 675 - 360(1 - Z2) + 16(1 - Z^f 



■) 1/2 

+ ((5 + 4(1 -Z2))(45 + 4(1- Z2))=')'/H . (27) 



As described above, we can obtain the 'absolute' value of ^Imax by setting Z = 0. This is 
the 'absolute' upper bound on the angular velocity which turns out to be 

Q'(l) - (2R-) 

A numerical value can be found by taking M — 1.35AfQ, so that 

f^:„,, « 2.35 X lOSad/s (29) 

is only one order of magnitude away from the observed millisecond pulsars. When Z is 
different from zero, ^max can be consistent with the largest angular velocities observed 
in millisecond pulsars The result for typical rcdshifts around Z ~ 1 is 

~ 2 ^ 



K..xiZ)^^[-^il-Z)) (30) 

One can interpret the equation ([50]) in two different, but related ways. Both ways have to 
do with the evidence that milli-second pulsars are orthogonal rotators {9 — tt/2) p fTUlITT] . 
Since the right hand side of (150]) agrees already with the angular velocity of fast spinning 
objects, the emission angle 9 must be close to tt/2 (see equation (IT4)) ). As we shall see 
below, this result can in turn be used to learn about the orientation of the magnetic axis 
in rapidly rotating objects, particularly neutron stars. 



The standard model for the pulsar emission mechanism was developed independently 
by Pacini [36] and Gold [37] (see also 'lighthouse' model [38]), and will now be described 
briefly. Since the rotation axis is not aligned with the magnetic axis, a changing magnetic 
field will induce electric flelds at the magnetic poles in for example, a neutron star. These 
electric fields will eject particles which will follow helicoidal paths around the magnetic 
field lines. The ejected particles will in turn, emit a narrow cone (~ 10° p[^) of radiation 
parallel to the magnetic axis. From this argument it can be inferred that the angle between 
the magnetic axis and the rotation axis is approximately the same as the emission angle 9. 
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Figure 3: Maximum angular velocity as a function of the redshift for different masses. 
For each mass, there are two curves. The upper curve is the maximum angular velocity 
(Equation (|27p ) obtained from the perturbative metric. The lower curve for each mass 
corresponds to the analogous expression for the extended metric (|33p . 



The emission angle has been measured indirectly for many pulsars by Kuzmin and Wu 
[5], and in their results, a strong correlation between orthogonal magnetic axes and fast 
millisecond pulsars is evident. This measured correlation agrees with our results and the 
results of [101 E], which favor orthogonal rotators. 

On the other hand, we can assume the theory outlined above to be valid, which allows 
us to determine ^Imax- Indeed, we can even assume that the angular velocity of fast ro- 
tating objects is close to the maximally allowed value. In other words we have f2 ~ ^max 
which essentially, knowing predicts the redshift. In turn, we can now extract the value 
of the radius. We will discuss this procedure taking realistic examples in the next section. 
Some examples of the behavior of S^'^ax ^ function of the gravitational redshift are 
shown in figure 3. 

There has been a considerable amount of work done related to the maximum angular 
velocity. Haensel et al. [13] , Glendenning 39J and Koranda et al. [3D] have used empirical 
formulae together with extensive numerical calculations to provide a lower bound for the 
period of rotation. Haensel et al. [13] found the minimum period to be Tmin = 0.288 ms, 
assuming a mass of 1.44M0. This period also depends on the minimum bound for the 
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redshift which was found to he Z = 0.528 [13]. If we use eq. (P7| (in accordance with 
the evidence of milh-second pulsars being orthogonal rotators we use 9 — 7r/2) and the 
parameters just mentioned we obtain Tmin = 0.442 ms. However, if we use our absolute 
bound (eq. and the mass assumed by Haensel we obtain the very close value of 

Tmin = 0.284 ms. This is just a difference of just 1.4% which is rather gratifying recalling 
that our result has been obtained analytically by methods different from [T5j . 

Other relevant work has been done related to the "mass-shedding limit" f[12j .[411 .[4211. 
This limit is conceptually different from what we have found and corresponds to the limit 
at which the neutron star would break apart. Lattimer and Prakash find for this limit 

^inaxiR) = 1045 {M lMQf/'^{lOkm/ Rf'^ Hz, (31) 

which is independent of the equation of state |17j and applicable for masses not very close 
the maximum mass. Equation (|3ip can be used for example, to find a limiting radius 
having an assumed or measured mass. While applying our results to find a radius we 
will compare our findings to the results when our maximal angular velocity is replaced 
by the one above. Here we note that that it is not straightforward to compare the radius 
independent limit (pS)) with pTjl as the latter depends explicitly on R. However, it is 
obvious that both are of the same order of magnitude. 

The orders of magnitude reached by the maximum angular velocity are quite large, but 
still less than the condition for slow rotation given by Zeldovich and Novikov [27]; even 
for the absolute maximum (Equation (|28p ). The condition for slow rotation (J <C MRgc) 
gives 

n « n^^L ^ (32) 

which is one order of magnitude larger than our maximum upper bound when we take R = 
Rg, and even larger when we compare with the non-zero redshift case. Still, because our 
maximum angular velocity is comparable to the Zeldovich-Novikov condition, especially 
with R > Rg, we will make a better approximation for the metric in the following section. 



5 Improvements: Extended Metric 

Even though the perturbative metric pop gives relevant results, a better approximation 
for the metric is given for instance in |43] . Assuming a negligible quadrupole moment, 
the second approximation the metric is 



9u 



9t4> 



J 2GM 1(GmV 
\ c^R 6\ c^R J 

Js\T?e i 2GM fGM\' 



GM 



?2 „:„2^ 



c^R 



\c^R ) 



(GM\ 



R:^sin\d) ] 2GM 1 fGMV 



O 



GM 



GM 



c^R 

3~1 



(33) 

(34) 
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One can apply the very same procedure as before. The resulting polynomial to be solved 
to obtain the two solutions for the radius is of fifth order in R. In consequence we 
can present the solution only numerically. An inspection of figures 2 and 3 shown that 
the new results, differ from the previous one only at very high angular velocities. The 
maximum angular velocity obtained using the extended metric is slightly higher than the 
one obtained analytically, so that eq. (177)) for the maximum angular velocity is still a 
valid bound. However, in comparison to the first order metric we used before, it becomes 
clear that the limit on the maximal angular velocity improves, especially for the absolute 
limit at zero gravitational redshift. 

6 Applications 

Here we shall apply our results to some known compact objects such as Sirius B, an iso- 
lated neutron star, and several known pulsars. 

Sirius B is the binary companion of the very bright Sirius A and is the closest White 
Dwarf to earth. Sirius B has been studied extensively, and the Gravitational Redshift has 
been measured accurately with the help of the Hubble space telescope [44]. Also, since it 
is a binary system, its mass has been measured accurately [45]. 

2 = 0.999735 ± 0.000015 ; M = 0.984M (35) 

According to eq. (j27p . we get a minimum period of 

iriusB) 



sm( 



= 10.8 s (36) 



This number does not change very much, had we applied the results from the extended 
metric. 



The low mass X-Ray binary system EXO 0748-676 has been studied by Cottam et al. 
[46] . They managed to measure spectral absorption lines corresponding to a redshift of 

Z = 0.74 (37) 

This together with its mass (Assumed to be 1.35M0) gives 

%n,n{EXOQ7AS ~ 676) ^ 9.05 x lO^^s (38) 

The value above has been obtained employing the extended metric (the corresponding 
value resulting from the first order approximation is 5.8 x 10~^s). The period for this 
neutron star has been measured to be 22 x 10~'^s[17], which is not too far away from 
the above limit given by Tmin ■ The predicted radius for the assumed mass and measured 
period and redshift, is 8.7km. 

In addition to the previous applications, our results concerning the maximum angular 
velocity, can be applied to millisecond pulsars. It is plausible to assume that for millisec- 
ond pulsars, the measured angular velocity is very close to the maximum angular velocity. 
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Pulsar 


Period (ms) 


M (M0) 


9 


Z 


R (km) 


PSR J1748-2446ad 


1.396 34J 


1.35? 


90? 


0.834 


20.1 


PSR B1937+21 


1.557 35J 


1.35? 


90 ^ 


0.836 


20.2 


PSR J1909-3744 


2.95 [48 


1.438 [48] 


90? 


0.896 


31.1 


PSR 1855+09 


5.3 


1.35? 


90 [g 


0.935 


46.9 


PSR J0737-3039 A 


22.7 [49 


1.34 [ig 


90? 


0.976 


133.6 


PSR 0531+21 


33.3 


1.35? 


90 9^ 


0.982 


164.4 


PSR B1534+12 


37.9 [W 


1.34 [Ml 


90? 


0.983 


165.9 



Figure 4: Predicted Maximum Z and radius for several pulsars. 



With this in mind, it is possible to predict a value for Z using equation (|30p . and a unique 
value for the radius using (fTT]) or ^8]) . Since Vlmax{Z) is a decreasing function, the pre- 
dicted Z is actually a maximum bound, and if the "preferred" radius for the neutron star 
is r2, the predicted radius can be thought of as a maximum bound also. Here we present 
a table with predicted redshifts and radii for several fast millisecond pulsars. 

From the table in can be inferred that for fast millisecond pulsars the radii are consis- 
tent with standard neutron star models |,19j . However, for the slower millisecond pulsars, 
the maximum radius is slightly greater than what is predicted by most neutron star mod- 
els, which implies that these neutron stars are probably not rotating exactly at their 
maximum angular velocity. In such a case the given radii should be interpreted as an 
upper bound which comes indeed close to the upper limit discussed in the introduction. 
Using equation ([3T|) and assuming a mass of 1.4 M© the corresponding radius for PSR 
B1937+21 turns out to be 15.5 km [17], 5 km less than our result. For the slower millisec- 
ond pulsars, our results also predict slightly larger radii than those obtained by equation 

7 Conclusions 

We have shown that the gravitational redshift in conjunction with global results from 
theoretical models can yield valuable information on the properties of rotating compact 
objects. Even though the determination of the radius in the presence of rotation is not a 
unique prescription, for relatively small angular velocities we can always opt for the lower 
result of the radius determination. With increasing angular velocity, a narrow range of 
possible radii can be defined. Alternatively, assuming that the angular velocity of the fast 
spinning neutron stars is close to its maximal value, we can either obtain a unique radius 
or an upper bound. The maximal angular velocity derived in text does not depend on 
the radius directly, but on the redshift which makes direct contact with existing or future 
observations. The absolute limit on Q. (|^^ does not even depend on the redshift. It is 
satisfying that both these bounds come close to the observed values for millisecond pulsar. 
This implies that nature reaches here it maximally possible value. Another advantage of 
our approach is the confirmation of the emission angle of radiation in fast rotating neutron 
stars. The application of our results to existing objects clearly show that the method of 
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using the gravitational rcdshift for rotating objects is effective. 

The important feature we would like to emphasize here is that we relate properties of 
Rotating Compact Objects to measurable quantities such as the Gravitational Redshift 
and the angular velocity. This way, our approach is semi-empirical and independent of 
model details. 

As discussed in section four, our analytical findings regarding the maximal angular 
velocity agree with results obtained after numerical calculations. In this way, both, the 
analytical and numerical approach, corroborate each other. 

Acknowledgment: We would like to thank Paolo Gondolo, from the University of 
Utah, for useful discussions and carefully reading the manuscript. 
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